Toward a microscopic description of dimer adsorbates on metallic surfaces 
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Despite the experimental successes of Scanning Tunneling Microscopy (STM) and the interest 
in more complex magnetic nanostructures, our present understanding and theoretical description 
of STM spectra of magnetic adatoms is mainly phenomenological and most often ignores quantum 
many-body effects. Here, we propose a theory which includes a microscopic description of the wave 
functions of the substrate and magnetic adatoms together with the quantum many-body effects. 
To test our theory, we have computed the STM spectra of magnetic Cobalt monomers and dimers 
adsorbed on a metallic Copper surface and successfully compared our results to recent available 
experimental data. 

PACS numbers: 72.15.Qm, 68.37.Ef, 72.10.Fk,75.30.Hx 



Introduction In the past decade, STM experiments of 
magnetic adsorbates on metal surfaces have florished. 
They have identified sharp resonances at the Fermi sur- 
face [E H S & @|, characteristic of the Kondo effect, a 
paradigmatic quantum many-body state which originates 
from the screening of the impurity spin by the conduction 
electrons. Systematic experimental studies of more com- 
plex structures with two or more magnetic atoms have 
been recently initiated and strong magnetic exchange in- 
teractions between the impurities have been identified 
d, 0] and estimated 0, 0| . Despite all these important 
experimental successes and the interest in more com- 
plex magnetic nanostructures, our present understand- 
ing is mainly phenomenological. Since the exponentially 
small Kondo energy scale cannot be captured by present- 
day calculations based on the Density Functional The- 
ory (DFT), a full ab initio theory of Cobalt impurities 
on metallic surfaces is precluded. A different strategy 
is therefore required. Here, we proceed in two separate 
steps by combining two complementary approaches: the 
single particle physics which is mainly responsible for the 
shape of the STM spectra is described by solving the 
Schroedinger equation in an effective potential which pa- 
rameterizes the metallic surface and the adatoms while 
the quantum many-body physics is treated separately 
via the powerful numerical renormalization group (NRG) 
method. Such strategy captures low energy physics and 
should therefore be suitable to describe STM spectra of 
magnetic adatoms around zero bias. This constitutes a 
first step toward a microscopic description of competing 
quantum many-body states at metallic surfaces. 

The Kondo effect occurs when the spin of magnetic 
adsorbate couples to the surrounding electrons. Be- 
low a characteristic energy scale Tk, named the Kondo 
temperature, the conduction electrons form a strongly- 



correlated singlet with the spin of the adsorbate which 
results in its screening. In STM experiments, this mani- 
fests as a zero-bias resonance of width Tk in the differen- 
tial conductance [E S & Hi ■ The theoretical description 
of the STM spectrum of a single magnetic atom adsorbed 
on a metallic surface is based on a Fano line shape as- 
sociated with a non-interacting resonant level. The Fano 
function is most often empirically assumed while its pa- 
rameters are adjusted phenomenologically 11, 12, l^, 14 1 
or computed from a microscopic description of mctal- 
adsorbate-tip interactions 1^, T^. 

When two magnetic adsorbates or more are brought 
into proximity, magnetic interactions between the impu- 
rities -direct or indirect- start competing with the in- 
dependent screening of each impurity. From this com- 
petition different magnetic ground states can be reached 
17, 3, lif that can strongly affect the differential con- 



ductance [20| and may even lead to non-Fermi liquid be- 
havior in small clusters [2l| . In Ref. , Wahl et al. per- 
form a systematic experimental analysis of Cobalt dimers 
adsorbed on a Copper surface by increasing the inter- 
atomic distance between the adatoms and therefore de- 
creasing their magnetic exchange interactions. The mea- 
sured changes of the STM spectra with the adatom in- 
teratomic distance thus permitted an indirect estimate 
of this coupling via STM. Below, we propose a theoreti- 
cal description of low energy STM spectra of magnetic 
monomers and dimers adsorbed on a metallic surface 
by combining microscopic evaluation of the one-electron 
parameters with NRG calculations without making any 
phenomenological assumptions on the spectral shapes. 
We focus on Cobalt atoms adsorbed on Copper (100) 
to compare our results with recent experimental data by 
Wahl et al. 5] but our approach can also be extended to 
other magnetic adatoms and metallic surfaces. 
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Model We consider a cluster of Nc magnetic atoms ad- 
sorbed on a metallic surface. The complete system in- 
volves the STM tip, the substrate and can be described 
by the following Hamiltonian 



H = Hs 



tip- 



(1) 



where Hgubs describes the substrate plus the adsorbed 3d 
transition metal atom and Hup-subs describes the inter- 
action of the tip with the substrate. Hup describes the 
tip which is assumed to have an unstructured density of 
states. The substrate with the adsorbed 3d atom may be 
modeled by a generalized Anderson model. The d^z^_r2 
orbital is the most strongly coupled to the metal at typi- 
cal adsorption distances of Co on Cu(lOO) due to the lobe 
pointing to the surface which overlaps more strongly with 
the metal wavef unctions. The d-orbitals parallel to the 
surface x-y (dxy, d^2_y2) plane lead to cancellations be- 
tween the negative and positive lobes of the adsorbatc 
orbital. This has been checked by evaluating the matrix 
elements along the lines shown in Ref. [15]. We therefore 
assume in the sequel that only the ^3^2^^^ orbital of the 
adatoms hybridizes with the substrate [22]. 

The metallic states, denoted by ]k >, couple to the 
1^322 orbitals denoted by ]di > with i = 1, .., Nc label- 
ing each adatom of the cluster. The Anderson model usu- 
ally assumes that the continuum of metal states are or- 
thogonal to the localized orbitals of the adsorbate. How- 
ever, the basis set formed by the unperturbed metal, ad- 
sorbates and tip states, {lk>,|c?i>,|i>}, is generically 
non-orthogonal and over-complete. One way to correct 
this is to redefine the metallic states, k, as: 
|k Ik > -X;^ < rfijk > \di > - < t\k > \t > . 
Considering that tip and adsorbate wavefunctions are or- 
thogonal: < t\di >— 0, as \di > is very localized, then 
the new metallic states satisfy < k]0 >= 0, where \4> > 
can be either \t > or \di >. 

Our starting point is thus an Anderson model defined 
in this new orthogonal basis: {]k >,\di >,\t >}, from 
which associated one-electron parameters are obtained. 
An analogous procedure was previously used in the con- 
text of chemisorption of atoms and molecules on metal 
surfaces by Grimley [2§]. In this basis, the substrate 
Hamiltonian reads 



Hsubs — e^'^t 



kcr jcr 



(2) 



j — l k.CT j,(7<(7' 

Here is the energy level of an electron residing in the 
d3,2_„2 orbital of the adsorbate, c\ creates an electron 

r _ ' kmc- 

with spin ct, momentum k in the metal, dj^ creates an 
electron with spin a in adsorbate-j. and Vj.^ are the 
metallic energies and the hybridization matrix elements 
between the substrate and the adsorbate-j, respectively 



(notice that the momentum dependence of the hybridiza- 
tion matrix elements is explicitly taken into account) . U 
is the Coulomb repulsion of two electrons in the 3d orbital 
of the transition metal atom. Finally, the tip-substrate 
interaction contribution to the Hamiltonian reads 



tip— subs 



(3) 



k,(T 



through the matrix elements, M^. Here, t„ destroys an 
electron with spin a in the tip. We have neglected the 
direct coupling of the tip with the substrate d bands and 
with the 3d orbital of the adsorbates due to the localized 
nature of these d orbitals. 

Hybridization matrix elements In the following we 
briefly describe how hybridization matrix elements, Mj^ 
and y^., are computed. For simplicity we first focus on 
how Vj. is computed, as Afj. is computed in a similar 
way. We first re-express the Hamiltonian in Eq. ^ in 

first quantized form Hgubs = T + X^^i ^ + ^M, where 
T is the kinetic energy of the system, Vj is the potential 
created by adsorbate-j and Vm describes the surface po- 
tential. The matrix elements between the orthogonalized 
metallic states, k, and the adsorbate thus read: 



Vkj - S'kj < dj\VM\dj >, 



(4) 



where, again, we have assumed: < t\dj >= and 
< d.i\dj >— Sij. The first term in Eq. ^ is the hybridiza- 
tion matrix clement with the unperturbed wavefunctions 
k: Vkj =< kjVA/lrfj >; ^-iid the second contains the over- 
lap matrix element S^j =< k]c?j >. The above orthogo- 
nalization procedure automatically selects the metal po- 
tential Vm in the hybridization matrix elements favoring 
the region close to or inside the metal in the integra- 
tions. This differs from hybridization matrix elements 
computed with the original wavefunctions as in that case 
integrations over the whole space are involved. The metal 
effective potential for the Cu(lOO) surface, Vm, is param- 
eterized based on the Jones- Jennings-Jepsen (JJJ) po- 
tential (see Ref. [13] for details). The sharpness of 
the surface barrier potential A = 2.2A~-'^, and the image 
plane position of Zim = 1-15 A. The adsorbtion distance 
of Co on Cu(lOO) is estimated to be Zd = 1.5 A. Finally 
we note that orthogonalization effects enter the model 
through matrix elements only. Orthogonalization effects 
on the substrate band energies can be shown to be of 
higher order in the overlap. Hence, we assume = ek 
in the rest of the paper. 



STM conductance results Following Ref. IJ] and for 
the sake of clarity we derive the basic equations needed 
for the computation of the conductance through the 
STM. If we neglect any modification of the substrate 
due to the presence of the tip (this is reasonable con- 
sidering the fact that the tip is typically at about 5 — 10 
A above the metal surface), then the differential conduc- 
tance measured by the STM reads: [III, [isl. [l4 [ 



3 



dV 



4e" 



-Ptip 



de^I^im+SVie)), (5) 



where r(a;) = TrJ^k. l-^kP'^('^~^k) is the conductance as- 
sociated with the clean substrate (without the adsorbed 
3d transition metal atom) and pup is the density of states 
of the tip and /(e) the Fermi-Dirac function. Modifica- 
tions of the tip-surface coupling induced by the presence 
of the adsorbate are given by 



i,j k,k' 



— Gij{uj) — 



IT] 



(6) 



The Greens function, Gij{uj), describes the electronic 
properties of the cluster of 3d adsorbates immersed 
in the metallic continuum including the quantum 
many-body effects such as the Kondo effect but also 
the interactions between the Co adatoms. rj is an 
analytical continuation parameter. For convenience 
Eq. ^ is rewritten in the following way 5r(w) — 
Im E,,, { (A, iLu)+iB, (w) )G,, (a;) (A* (u;) + iB* (u;)) } where 
Bi{uj) = vr J2k ^^i^k^i^ ~ ^k) and AiiiJ) is the Kramers- 
Kronig transformation of BiiiS). For the systems of in- 
terest here, Ai(uj) and Bi(io) are real. Function Bi(uj) 
embodies the information concerning the tip-substrate- 
adsorbate system as it depends on the tip-adsorbate sep- 
aration, on the position of the adsorbate with respect to 
the last surface plane of ions, and on the metal poten- 
tial described by Vm through the matrix elements, 
and Afj^. In the single impurity case, dI/dV{uj) reduces 
to a Fano form[l3| for uj < Tk and the line shape is 
determined by the Fano parameter [is}: q = A{Q)/ B{Q), 
evaluated at the Fermi energy taken as zero. Thus, the 
shape is fixed by the one-electron parameters entering 
the full model 

We apply this formalism to compute the conductance 
of two Co atoms deposited on Cu(lOO) surfaces. We take 
in the following the origin of distances at one of the Co 
adatoms in the dimer Ri ~ and the distance between 
the atoms is denoted by d. The magnetic exchange in- 
teraction depends on the interatomic distance d and is 
denoted as I{d). This function has been calculated using 
ah initio calculations by Stepanyuk et al.[25| and also 
estimated experimentally Due to the small spatial 
range of the d orbitals compared to the typical separa- 
tion between the adatoms considered here, this magnetic 
exchange interaction is mainly of RKKY type as it has 
been shown in (25l |. For our purpose, we merely regard 
I{d) as an external parameter which may eventually ex- 
tracted from ab initio calculations. 

The challenging task is to compute the two by two ma- 
trix Gij in Eq. ([6|) which encodes interaction effects. In 
order to perform this task, we use the powerful NRG ap- 
proach. However, a full NRG treatment of two five-fold 



degenerate impurities is a way beyond today's computa- 
tional possibilities. We choose instead a simpler route 
and first map the Anderson Hamiltonian in Eq. ([3]) to 
a Kondo-like Hamiltonian. Such procedure is standard 
and justified by the fact that the electron in the d^z^_r2 
orbital of the Cobalt adatom is strongly bound. The re- 
sulting Kondo model is described by the two-impurity 
Kondo Hamiltonian [l3, • We neglect potential scat- 
tering terms since they do not renormalize and barely 
contribute to the low energy physics. The low energy 
physics of the two-impurity Kondo model is then con- 
trolled by only two energy scales: the single-impurity 
Kondo tem pera ture Tr- and the magnetic exchange in- 
teraction / 17, ISi. We use NRG calculations to obtain 



the impurity spectral properties. 

Wilson's NRG technique [i^ is a non-perturbative and 
numerically exact method suitable especially for quan- 
tum impurity problems (for a recent review see [13 )■ The 
cornerstone of the method is the logarithmic discretiza- 
tion of the conduction electron band and mapping the 
system onto a semi-infinite chain with the impurity at the 
end. NRG has been successfully applied to the problem 
of two magnetic impurities [l^, [l^ but in those earlier 
works only thermodynamic quantities were computed. 
The main source of complication in a two impurity cal- 
culation as compared to a single impurity case is that the 
former is an effectively two band calculation: From NRG 
point of view it is a challenging task because the impuri- 
ties now couple to two semi-infinite chains. Consequently, 
the Hilbert space grows by a factor 16 in each NRG step. 
This is still manageable with today's computer resources. 
Concerning the details of the present NRG calculation, 
the conduction and is discretized logarithmically in in- 
tervals [A"+ii:), A"D] and [-A"D, -A^+^D] where D is 
the bandwidth and A the discretization parameter. We 
took A = 2, the number of iterations was A'' = 60 and we 
kept M = 3072 states per iteration exploiting the charge 
and spin z-component U{1) symmetries. The calcula- 
tions were performed at T = and /(rf) was taken as an 
input parameter we took from ah initio calculations [25| . 
In the NRG calculations, we have neglected the details of 
the substrate band structure and the ^-dependence of the 
matrix elements Vkj since the low energy properties of the 
two impurity Kondo model mainly depend on the single 
impurity Kondo temperature and the RKKY interaction. 
Therefore, we expect our whole approach combining sin- 
gle particle description of the adsorbate together with 
NRG calculations to be qualitatively correct at low ener- 
gies i.e. around zero bias and within an interval of order 
a few Tk. It would be interesting but quite computa- 
tionally demanding to apply NRG techniques directly on 
the two-impurity Anderson Hamiltonian neglecting or- 
bital degeneracy with the matrix elements Vkj computed 
from the previous single-particle calculations. We leave 
this for future work. 

We would like to note here that in principle NRG is 
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FIG. 1: Kondo resonance for a single Co adatom on a Cu(lOO) 
surface. STM measurements at T=6 K (red line) are com- 
pared to our theoretical results (black line). 



capable of reproducing RKKY interaction, see Ref. [19|. 
Here we choose a different strategy, by mapping the prob- 
lem onto two independent Wilson chains, but in contrast 
to Ref . [ist , we take the finite overlap between the elec- 
tron fields centered around the positions of the two im- 
purities into account. In this way NRG does not generate 
additional RKKY interaction thus the danger of double- 
counting is avoided. The reason why we follow this strat- 
egy is that the RKKY interaction, generated in the NRG 
scheme of Ref. [19] for a Kondo-like model with constant, 
featureless conduction electron DoS, has a very little to 
do with the realistic 1(d) obtained by ab initio method. 

When the two Cobalt adatoms are far apart, the 
RKKY interaction is negligible and the adatoms can be 
regarded as isolated. We have plotted in Fig. [T] the dif- 
ferential conductance -function of the dimensionless pa- 
rameter Lo /Tk- calculated using the aforementioned pro- 
cedure for a single Co adatom and compared it to exper- 
imental raw data. The agreement is extremely good at 
low energy and qualitatively correct at higher energy. 

Let us now apply the same method to obtain STM 
spectra for Cobalt dimers separated by a distance d such 
that the RKKY interaction become of the order of Tk 
or larger. The resulting differential conductance on top 
of one adatom are summarized in Figs. [2] and [3] for var- 
ious distances separating the Cobalt adatoms. Accord- 
ing to ah initio results [25|, we can associate to the dis- 
tance A—b.72 A and d=5.12 A the antiferromagnetic 
RKKY exchange interactions / ~ Tk « 7.6 meV and 
/ ~ 2Tk ~ 15.6 meV respectively. We have compared 
in Fig. [2] our results for the conductance on top of one 
adatom for these two distances to the experimental raw 
data obtained at T = QK with a very good agreement 
especially at low energies. The agreement is quite re- 
markable since the theory contains a single fit parameter: 
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FIG. 2: Kondo resonance for Co antiferromagnetically cou- 
pled Co dimers on a Cu(lOO) surface. The distance, d, (in 
A) between the Co atoms is reduced from (a) to (b) with the 
corresponding coupling /(d) [2^ increasing as shown. In (a) 
and (b) we compare theoretical calculations for an antiferro- 
magnetic interaction, / > 0, and available experimental data 
at 6 K. In (b) a linear background present in the experi- 
mental data has been appropriately added to the theoretical 
spectra for comparison. 



the magnetic exchange interaction / which is extracted 
from ab initio calculations [25|). The correct shape of the 
STM signal is reproduced without using any Fano line 
shape fits which would be actually unjustified for Cobalt 
dimers. We then apply our method in order to predict the 
differential conductance for the adatom separation cor- 
responding to d=3.5 A and d=2.56 A as shown in Fig. 
[3l The magnetic exchange interaction is predicted to be 
ferromagnetic [2^ and corresponds to / w —4Tk) and 
/ w —AQTk respectively. In the latter case, the Kondo 
temperature is predicted to be very small, « 0.23 meV 
and although is not accessible with present STM exper- 
iments it is the predicted line shapes resulting from our 
combined approach at very low temperatures. 

In summary, we have presented theoretical predictions 
for the STM spectra of magnetic monomers and dimers 
adsorbed on metals. Our approach relies on a micro- 
scopic description of the adsorbate-metal-tip interactions 
combined with NRG calculations. Such methodology 
may constitute the first building block for understand- 
ing STM spectra of more complex magnetic structures 
on metal surfaces. 
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FIG. 3: Kondo resonance for ferromagnetically coupled 
dimers on a Cu(lOO) surface. Experimental data do not show 
a Kondo effect at 6 K for d=2.56 A, which is in agreement 
with the strong reduction in the Kondo temperature predicted 
here. 



[1] H. C. Manoharan, C. P. Lutz, and D. M. Eigler, Nature 

(London) 403, 512 (2000). 
[2] W. Chen, T. Jamneala, V. Madhavan, and M. F. Crom- 
mie, Phys. Rev. B 60, 8529 (1999). 

T. Jamneala, V. Madhavan, and M. F. Crommie, Phys. 
Rev. Lett. 87, 256804 (2001). 

C. F. Hirjibehedin, C. P. Lutz, and A. J. Heinrich, Sci- 
ence 312, 1021 (2006). 
[5] P. Wahl et ai, Phys. Rev. Lett. 98, 056601 (2007). 
[6] V. Madhavan, W. Chen, T. Jamneala, M. F. Crommie, 
and N. S. Wingreen, Science 280, 567 (1998). 
J. Li, W.-D. Schneider, R. Berndt, and B. Delley, Phys. 
Rev. Lett. 80, 2893 (1998). 
[8] N. Knorr, M. A. Schneider, L. Diekhoner, P. Wahl, and 

K. Kern, Phys. Rev. Lett. 88, 096804 (2002). 
[9] N. Neel, J. Kroger, L. Limot, K. Palotas, W. A. Hofer, 



[3] 
[4] 



[7] 



[10 
[11 

[12; 

[13 
[14 

[is: 

[16 
[17 

[is: 

[19 

[2o: 

[21 

[22: 
[23: 

[24: 

[25: 

[26' 
[27 



and R. Berndt, Phys. Rev. Lett. 98, 016801 (2007) 
U. Fano, Phys. Rev. 124, 1866 (1961). 

A. Schiller and S. Hershfield, Phys. Rev. B 61, 9036 

(2000) . 

O. Ujsaghy, J. Kroha, L. Szunyogh, and A. Zawadowski, 
Phys. Rev. Lett. 85 2557 (2000). 

M. Plihal and J. Gadzuk, Phys. Rev. B 63, 085404 

(2001) . 

V. Madhavan, W. Chen, T. Jamneala, M. F. Crommie, 
and N. S. Wingreen, Phys. Rev. B 64 165412 (2001). 
J. Merino and O. Gunnarsson, Phys. Rev. B 69, 115404 

(2004) ; Phys. Rev. Lett. 93 156601 (2004). 

C.-Y. Lin, A. H. Castro Neto, and B. A. Jones, Phys. 
Rev. Lett. 97, 156102, (2006); Phys. Rev. B 71, 035417 

(2005) . 

B. A. Jones, CM. Varma, Phys. Rev. Lett. 58, 843 

(1987) . 

B.A. Jones, CM. Varma, and J.W. Wilkins ibid. 61, 125 

(1988) . 

J. B. Silva, W. L. C Lima, W. C. Oliveira, J. L. N. Mello, 
L. N. Oliveira, and J. W. Wilkins Phys. Rev. Lett. 76, 
275 (1996) 

P. Simon, R. Lopez, and Y. Oreg, Phys. Rev. Lett. 94 
086602 (2005). 

B. Lazarovits, P. Simon, G. Zarand, and L. Szunyogh, 
Phys. Rev. Lett. 95, 077202 (2005); K Ingersent, A. 
W. Ludwig, and I. Affleck, Phys. Rev. Lett. 95, 257204 
(2005). 

M. Weissmann and A. M. Llois, Phys. Rev. B 63, 113402 
(2001). 

T. B. Grimley, Molecular processes on Solid Surfaces, 
edited by E. Drauglis, R. D. Gretz, and R. L Jaffee 
(McGraw-HiU, New York, 1969), p. 299. 
R. O. Jones, P. J. Jennings, and O. Jepsen, Phys. Rev. 
B 29, 6474 (1984). 

V. Stepanyuk, A. Baranov, D. Bazhanov, W. Hergert, A. 

Kastelson, Surface Science 482-485, 1045 (2001). 

K. G. Wilson, Rev. Mod. Phys. 47, 773 (1975). 

R. Bulla, T. A. Costi, T. Pruschke, Rev. Mod. Phys. 80, 

395 (2008). 



